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Abstract 

This paper concerns the reconstruction of the scattering coefficient in a two-dimensional 
transport equation from angularly averaged measurements when the probing source is isotropic 
and time-harmonic. This is a practical setting in the medical imaging modality called Optical 
Tomography. As the modulation frequency of the source increases, we show that the recon- 
struction of the scattering coefficient improves. More precisely, as the frequency w increases, 
we show that all frequencies of the scattering coefficient lower than b are reconstructed sta- 
bly with an accuracy that improves as lo increases and b decreases. The proofs are based 
on an analysis of the single scattering singularities of the transport equation and on careful 
analyses of oscillatory integrals by stationary phase arguments. 

1 Introduction 

In the theory of stationary inverse transport, it was shown in [7] that, with knowledge of the 
attenuation coefficient, the problem of reconstructing the spatial part of the scattering coeffi- 
cient k{x) from isotropic sources and angularly averaged measurements was severely ill-posed. 
This was due to the fact that the measurements were not capturing enough singularities of the 
transport solution because of angular averaging. In the present paper, we consider a more favor- 
able measurement setting where the (still isotropic) source term is modulated in frequency with 
a sufficiently high modulation frequency. For sufficiently high frequencies, this corresponds to 
measurements of the time-dependent transport equation at very short time sampling (theory for 
this can be found in [5]). In either case, asymptotic expansions (for large modulation frequencies 
or small times) of the single scattering component of the measurements show that in the leading 
order, and in the case where the domain of interest is a ball, we recover a weighted X-Ray 
transform of the scattering coefficient O [6] , while the remainder term in the measurements is 
of lower-order. The practical interest of using time-harmonic sources in transport for improving 
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the reconstruction of optical parameters was first pointed out numerically in [T3] in the context 
of least-squares based reconstructions, where the authors observed that cross-talk between re- 
constructed optical coefficients was reduced and accuracy of the reconstructions was improved 
as uj increased. One of the goals of the present paper is to provide theoretical justifications for 
these observations. 

Although the remainder in the asymptotic expansion of forward transport solutions has small 
magnitude (see [6]), it can still be magnified during the inversion procedure because the inverse 
X-Ray transform is a deregularizing operator. In order to avoid this, we need to regularize the 
inversion procedure. We do so by giving up on the reconstruction of the high-frequency content 
of the unknown scattering coefficient k. Our main objective is to show which spatial frequencies 
of k may be stably reconstructed and with which accuracy as a function of the modulation 
frequency u. 

The task of the present paper is to devise a proper reconstruction algorithm for the scattering 
coefficient in this setting that is justified theoretically. The theory is based on stationary phase 
arguments, which allow us to exhibit decay of solutions (see also [6]) and errors in reconstructions 
in terms of the modulation frequency cj. The same techniques are also used to guarantee the 
convergence of an iterative reconstruction scheme and in particular to prove that certain error 
operators become contractions when oj becomes large enough. 

The results presented here show that in order to attain a given precision on the reconstruc- 
tion, there is a trade-off between the modulation frequency uj of the measurements available, 
and the bandwidth b (or fineness of details) at which one will recover the function k. Knowledge 
of measurements at a given frequency uj allow for a stable reconstruction of the scattering coef- 
ficient up to a maximal resolution b that depends on this uj and the accuracy we wish to attain. 
Modulated sources at a given frequency corresponds to a spatial scale below which details of k 
are not accessible. We make such statements precise. 

The rest of the paper is structured as follows. We state our main results in section [2j Section 
13.11 recalls the necessary elements of forward transport theory in a general setting, including 
the asymptotic behavior of angularly averaged solutions of time-harmonic transport for high 
modulation frequencies. Section 13.31 describes how to invert for the spatial part of scattering 
k{x) in two dimensions of space and when the domain is a disk. Sections 13.41 and H] provide 
proofs for the main results. 



2 Statement of the main results 

Let r > and let Br denote the open disk {xf + < r^} in M? with boundary denoted by dBr- 
We consider the two-dimensional time-harmonic transport equation on B,. x 

J§i (2.1) 

u{x, v) = 5{x — xq), {x, v) G r_, 
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where 6{x — xq) is an isotropic point source located at xq and is defined for any test function tp 
on dBr by Jq^ 6{x — xo)ip{x) dfi{x) = ip{xo). 

In this paper, we consider angularly averaged outgoing measurements given by 

T^{xo,Xc) ■= / u\Y^{xc,v) V ■ dv, u solves ()2.ip . (2-2) 

Such measurements are reasonable descriptions of the data acquired in practice in the medical 
imaging modality called Optical Tomography; see [13] and references there. Because the source 
term in (|2.ip is isotropic and the measurements are angularly averaged, the main singularities 
of the transport equation are not captured by such measurements. The reconstruction of the 
optical coefficients thus becomes a severely ill-posed problem [H [7]. One way to restore some 
well-posedness in the inversion is to consider the limit as w — )• oo. This is the problem considered 
in this paper following results obtained in [6] on the forward transport problem. 

Under appropriate subcriticality conditions recalled in the next section, we obtain existence 
and uniqueness of the solution of (j2.ip as well as the decomposition 

u{x, v) = J^5x^^ + IC^JuiS^o + {1 - ^uy^K-l^Juj^xo, (2.3) 

where the operators J^) and K.^) are defined in ()3.6p and (j3.5p below, respectively. This decompo- 
sition in turn yields a decomposition of the measurement operator (12. 2p into three components 
denoted by 

[a, k]{xo,x,):= [a] {xq , x^) + [a, k] {xq , x^) + T!^+ [a. A;] (xq , x,) , (2.4) 

where Tq accounts for particles that traveled from xq to Xc in a ballistic way and Ti and T2+ 
account for particles that were emitted at xq, scattered once or multiple-times inside the domain, 
respectively, and were measured at Xc- Here and below, xq denotes the emitter's position and 
Xc, the detector's position. 

In this paper, as in [7j, we assume a to be known and focus on the reconstruction of the 
scattering coefficient k{x). The ballistic part can be taken out of the measurements and we 
define our data to be 

D'^[/c](xo,Xc) ■.= T^[a,k]{xo,Xc)-T^[(j]{xQ,Xc) = T^[k]{xo,Xc) + T!^+[k]{xo,Xc). (2.5) 

In the sequel, we address the reconstruction of the function k from knowledge of the data ()2.5p . 

Throughout the paper, we require the following crucial hypothesis, which expresses the fact 
that k vanishes in the vicinity of dX: 

Hypothesis 2.1. There exists < D < r such that supp k C 5,-21) • 

The reason is that the measurements we consider are overwhelmed by a term that depends 
only on the value of k at the domain's boundary. In other words, to leading order, the mea- 
surements do not depend on the scattering coefficient k away from the boundary. We thus 
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assume here hypothesis I2.H which in practice corresponds to placing the detectors away from 
the scattering domain of interest and thus does not seem to be overly restrictive. 

Now, as is proved in [6] assuming a hypothesis equivalent to 12.11 the expression ()2.5p admits 
the following asymptotic decomposition for large w: 



D'^{xo,Xc) = Tf'[k]{xo,Xc) + Tiji[k]{xo,Xc) + T2j^[k]{xo,Xc), where 

1 f C (2-6) 



— If c 

T^[k]{xo,Xc) = -^A^{xo,Xc) / kp, and ||T{^r[A;] + T2+[k]\\ ugx)2\ < — , 



and where is bounded away from zero independently of uj and p{x) = (r^ — |a;p)~2 ^ x £ Br- 
In other words, the single scattering decomposes into a leading part that is proportional 
to the X-Ray transform of kp (and from which we would like to reconstruct all or part of k) 
and a remainder Tj^ that is asymptotically smaller. Equality (j2.6p tells us that for large uj, the 
remainder + becomes negligible. 

In order to setup an inversion, we first reparameterize the data into the parallel geometry 
of the X-Ray transform, namely we define Xq G dBr and G dB^ as functions of {s, 0) E 
(— r, r) X S-"^, such that the line joining xq and Xc is exactly 



L{s,e) := \^se^ + te, teM}, 



and rewritting equality (12. 6p in this parametrization yields 

x{sW{s,e) = ^A^{s,e)p[pk] {s,e) + -x{s)R^{s,e), (2.7) 



where P denotes the X-Ray transform in the parallel geometry, and x is a smooth function 
satisfying x{s) = 1 if |s| < r — 2D and xi^) = if |,s| > r — D. As will be seen later, x is 
necessary in order to smooth out effects at the boundary of the term , while leaving the term 
P [pk] unchanged thanks to hypothesis 12.11 

Let us now set up an inverse for T|^[fc]. The first candidate for an inversion is 

Tf^' = ^P-^A')-^.], 

where P~^ denotes the inverse X-Ray transform. Such an inverse is unsuitable because P~^ is 
a deregularizing operator that prevents us from controlling the error term o [T{^ -|- T^] 

in a suitable manner. We therefore introduce a regularized inverse X-Ray transform P~^'^ with 
6 > 0, such that 

p-^'''oP[f]{x)=W,^f{x):=h{x), where W,{x) := ^ j^^e'^'^^ {^^^ (2.8) 
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with ^> : [0,00) I— 7- [0,1] a "low-pass filter" supported inside [0,1]. Hence, the parameter ^ 
measures the finest scale at which we reconstruct / when applying P~^'^ to P[f]. The inverse 
operator that we apply to the data (j2.5p is therefore defined by 



Tf := V^p-i'^[(A-)-ix-]. (2.9) 
As the operator ()2.9p is applied to the data (I2.5p . we obtain an equation of the form 

f^~^'^T)^[k] = [kp]b + R'^'''[k], (2.10) 

where the "error" term R^'^lk] is decomposed into i?*^'* = R!^'^ + i^g'^'i with 

Rf[k] := Tr~''Vi?W = Tr~''V[^] - [kp]b, (2.11) 
i?2'''[A;] := 2f"^'V2"'+[A:]. (2.12) 

Here, R^'^lk] represents the error due to the part of the single scattering that is not inverted for 
whereas R^^'^lk] represents the error due to multiple scattering. 

Our first result controls the remainder R'^'^[k] in L°°{Br) using the asymptotic decomposition 
for large cj of the measurements (j2.6p as established in [6]: 

Theorem 2.2 (Direct inversion). Let a E C^{Br) be known and assume that k G C^{Bj.) and 
that the support hypothesis \2.1\ holds. Then for u large enough and fixed b, knowledge of the 
data T)^{s,6) for s G [— r + D,r — D] and 9 £ allows us to reconstruct [pk]ti up to an error 
i?^'*[A:] (fZTO]) that is bounded in L'^{Br) by 

\\R''''[k]\\L^<C-^. 



Here, the constant C depends on \\k\\Q2, \\a\\Q2, D and \\wi\\ii. 

Remark 2.3. Theorem \2.^ merely states that in order to get details of size ^ while keeping 
a constant error level, we must provide measurements at frequency uj ^ b'^ . Conversely, if 
measurements are available at a given frequency uj, then theorem \2.2\ gives the increase in error 
as we try to recover finer details of [pk] . 

In the second result, which is the main result of the paper, we assess the conditions under 
which the following iterative algorithm: 

ko = TY'"^'^D'^, 

^ (2.13) 



kn+i = G{kn) := ^'"2)'^ - R'^'^ikn], n > 
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improves the reconstruction of [kp\i, for a given b. The reconstruction in Theorem 12.21 treats 
multiple scattering as an unknown error term. However, once an approximation of k is obtained, 
we can estimate the multiple scattering contribution to the data and remove it from the data 
is order to obtain a better approximation of the single scattering contribution on which the 
reconstruction algorithm is based. Improvements are then obtained iteratively as described in 
(I2J3]) . 

Our objective is to show the improvements resulting from using such an iterative scheme 
and to obtain conditions that guaranty its convergence. In particular, we must find a suitable 
space for k and a (6, a;)-regime where the operator G is a ci-contraction for < ci < 1. This in 
turn requires a thorough and somewhat surprisingly complicated analysis of the operators R^' 
(j2.1ip and i?2 (j2.12|) that includes a careful analysis of oscillatory integrals by the method of 
stationary phase with errors that depend on the function k but not on its derivatives. 

The following theorem states that such a regime exists and details the accuracy that is 
reached in this case. In particular, we emphasize the fact that no regularity conditions are 
imposed on k. 

Theorem 2.4 (Iterative improvement). Assume that a G C'^{Br) is known and that k satisfies 
hypothesis \2.1\ Suppose further that the functionWi in (I2.8P satisfies ||Tyi||/^i(K2) < oo. Then for 
fixed b > 0, there exists Ki > and ojq > 1 such that for uj > coq and k such that \\kp\\oo < Ki, 
the iterative scheme (|2.13p converges to k* e BKo{L°°{Br-.2D))- Moreover k* satisfies the error 
estimate 

\\k^ - [kpUoc < T^\\[kp] - [Mfelloo, (2.14) 
1 - Cl 

where c\ G (0, 1). 

For fixed 6, the constant c\ behaves in uj like 

ci-^(6^ + 6^) + C,4log(^^), 

where the first and second terms in the right-hand side are constants of boundedness in L°° {Br-2D) 
for the error operators B!^' and R2' , respectively. 

Remark 2.5. Theorem \2.4\ implies a much better accuracy at the cost of a more expensive 
inversion because of the iterative scheme. At each iteration, the computation of R^'^ requires 
the computation of a forward transport problem. 

Remark 2.6 (On generalizations to higher dimensions). The extension of Theorem \2.S\ to higher 
dimensions should cause no difficulty as the proof of this theorem mostly relies on the asymptotic 
expansion of the forward solutions derived in for all dimensions. The extension of Theorem 
\2.4\should prove much more challenging. The fixed point method applied in Theorem\2.4\ requires 



6 



that the stationary phase integrals compensate the inversion of a term linear in k whose amplitude 

n — 1 

of order u 2 grows with spatial dimension n. In dimension n = 2, the stationary phase integrals 
provide 0{u!~^) contributions that compensate the previous term. The difficult problem of how 
the stationary phase techniques, which are at the core of our proof of Theorem \2.4[ should be 
modified in higher spatial dimensions is not considered here. 



Roadmap of proofs : Theorem 12.21 almost comes for free once the results from [6] are stated 

since it amounts to bounding the composition of bounded operators Tf o (T{^ + T^j. On 
the other hand, theorem 12.41 involves a much more technical analysis that will be the main focus 
of the paper. The proof relies on Propositions 14.21 and 14. 3^ which give bounds of the error 
operators R^'^ and R^'^ in some C[L^{Br^D)) spaces, respectively. For fixed values of 6, these 
norms decay with uj and so these operators become contractions for large enough u. 

Both propositions require non-standard estimates on oscillatory integrals (01) sometimes 
involving caustic sets (degenerate critical points of the phase function) near which one has to 
refine the analysis. References [U El [121 E] have greatly influenced the approach adopted here. 
Since proposition 14.31 involves technical analysis as well, we have postponed the stationary phase 
(SP) part of this proposition to proposition 14.51 Finally, in order to avoid redundancies in the 
paper, we have formulated an auxiliary lemma WA{ which gives, under a non-degeneracy condi- 
tion, an accurate first-order expansion of a parameter-dependent 01 with a simple expression of 
the remainder. Lemma l4.ll is used to prove propositions 14.21 and 14.51 



3 Forward transport theory, setup for inversion, and proof of 
theorem 12.21 

3.1 Forward transport theory 

We now recall some general results about forward transport that are necessarily in our analysis. 
Let X C M",n > 2 be an open convex bounded domain with boundary dX and diameter 
A > 0. Denote the incoming and outgoing boundaries 

r± = {{x,v) edX X S"-^ \±u^-v>0}, (3.1) 

where Vx is the outer normal to dX at x G dX. The time-harmonic transport equation reads 

V ■ Vu + {a{x,v) + iijj)u = / k{x,v' ,v)u{x,v') dv' , (x, f ) G X x S"^"*^, 

j§"~i (3.2) 

u{x,v) = g{x), {x,v) G r_, 

where w > and the input function g takes the form g{x) = S{x — xq), {xq,x) G (dX)'^ (call 
it )■ By 5{x — Xq) we mean the delta distribution that satisfies for each smooth function tp 
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defined at the boundary: 



5{x — xq)'iP{x) dfi{x) = ip{xo), (3.3) 

dx 

where d^{x) is the intrinsic measure at the boundary. For {x,v) G (X x S"^i)Ur+(jr-, 
let T±{x,v) be the distance from x to dX traveling in the direction of it-y, and x±{x,v) = 
X lb T±{x, v)v be the boundary point encountered when we travel from x in the direction of if. 

As it is done in many settings, we integrate (|3.2|) along the direction v. We obtain that u is 
a solution of the following integro-differential equation 

(I - IC^)u = J^g, (3.4) 

where we have defined, for G L^(X x S"^^) and ij) G L}{dX) 

IC^(p{x,v) := e~''^^E{x -tv,x) k{x - tv,v' ,v)(l){x - tv,v')dv' dt, (3.5) 

JO J§"-i 

J^^p{x,v) :=e-''^^-^'''''^E{x-{x,v),x)^p{x-{x,v)), and (3.6) 
E{x, x) := exp — / a{x + x' — xs, x' — x) ds \ . (3.7) 



The operators /Coj and J^j are well-defined and continuous operators in C{L^{X x ^)) and 
C{L^{dX),L^{X X §"-^)), respectively [HI [TO]. 

We now state a theorem recently proved in [6], which exhibits the asymptotic behavior for 
high oj of the solution u of (13. 2p . and motivates an inversion formula for k ([•] denotes the ceiling 
part of a real number): 

Theorem 3.1. Assume that {a, k) G C^^^ (X x S''-^) x C^^^ {X x §""1 x S"^i), and assume 
hypothesis \2.1\ Then, for large enough to, the measurement function T'^ admits the following 
singular decomposition 

T'^{xo,Xc) = T^ixQ, xc) + T^ixo, Xc) + T^+ixo, Xc), (3.8) 
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for a.e. {xo,Xc) € dX x dX, where 

Tg (xo,Xc) = ■. — - — r;;—^ \iy^o ■ (^olWxc -eol, (3.9) 

Ti{xo,Xc) = (xo, Xc) + T^r{xo, Xc), where 

( n — l ^ 
^) ' E{xo,x,)W,,-eo\Wx^-eo\ /" ° M^o + ^eq, eo, eg) (3^^^^ 



n+l 



Ti^R G L°^(5X X 9X) and ||ri"}j||oo < Ci (3.11) 

C2 <^~^, n = 2, 

Ta^^ G L°^(5X X aX) and \\T^+\\oo < { C2Uj-^Huj), n = 3, (3.12) 

(72 w 2 , n > 4, 



where Ci and C2 depend on the t ~^ and t ~^ norms of the coefficients a, k, respectively. 

In two dimensions of space, theorem 13.11 req uires three derivatives on the optical coefficients 
to ensure that decays hke , and only two derivatives to ensure that T^^ decays like uj~^ 
in L°°{dX X dX). Since the total remainder Tj^ + T^ cannot decay faster than w"^, we choose 
the a;~^-decay of and save one derivative on the optical coefficients. Theorem 13. II becomes: 

Theorem 3.2. Let n = 2, and assume that {a,k) G C'^(X x S"~i) x C^(X x S"~i x S""^) and 
that hypothesis \2. 1\ holds. Then the decomposition of theorem \3. 1\ holds with the following bound 

ll^lijlloo + ||^2+lloo ^ ^, (3.13) 

where the constant C depends on the norm of the coefficients (a, k) . 

Proof. In the proof of [U Lemma 5.2], do not perform the last integration by parts in the case 
m + d odd. It trades one derivative of the integrand for a decay factor of a;~2 . □ 

3.2 The 2D case, when X = Br and k(x,v,v') = k{x)(j){v,v') 

In this paragraph, we assume that the scattering function admits the expression k{x,v,v') = 
k{x)(j){v,v'), where </> is a known phase function. Suppose that the domain X is the centered 
ball of radius r > in M". For xq and Xc located on dX, we have the properties 

{xc - xq) ■ eo = \xc - xo\ and (xc + xq) • eo = (l^d^ - |xo|^)|xc - xoT"^ = 0, 

from which we deduce that — 2xo • eo = \xc — xq\. In this case, 



r — |xo + teo| =r — |xo| — 2txo • eo — t =t\xc — xo\ — t =|xc — xo|t 1 
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Hence we can rewrite the first term in equation p.lOp as 

Wou^J Jo - + teo\'^) — 

(3.14) 

In this case, the hne integral in ()3.14p is no longer a general weighted integral transform but 
simply the X-Ray transform of the function k{x){r'^ — |xp)^~, whose inversion is much more 
practical. Back to the two-dimensional case, we can now deduce the expression of the function 
A'^ defined loosely in 

A'^{xo,x,) = e~''^'^" (J^y Eixo,Xc)W^, ■ eo\W,^ ■ eo\Heo,eo)- (3.15) 

3.3 Inversion in two dimensions in parallel geometry 
3.3.1 The parallel geometry 

In order to setup an inversion, we must now choose a geometry that is adapted to the X-Ray 
transform. Namely, we must introduce variables that parameterize lines in space. We choose 
the parallel geometry since inversion formulas in this case are suitable for our analysis. 

In parallel geometry, we introduce the variables (s, 6) £ Z ■.= M. x (0, 2tt) and we would like 
to parameterize the emittor's position xq and the detector's position Xc in such a way that the 
line joining xq and Xc is 

L{s,e) := \^s9^ + t9, tGM}, 

where we have defined 6 = (cos 0, sin 0) and 6^ = {— sin , cos 9) . With the constraint that 
l^^ol = \xc\ = 1^, this is done by writing 

xo{s, 9) = s9^ - \/r2 - s'^9, and Xc{s, 9) = s9^ + \/r2 - ^2^, (3.16) 

in which case we have eo := Xc — xq = 9 and do{s) := \xc — xo\ = 2\/r'^ — s^. 

We now denote by = [—r + D,r — D] x the support of x{s)D'^{s,9). Over Zd, we 
have that Vr'^ — s"^ > V rD > 0. This guarantees that xq, Xc, do{s) and eo(^) are well-defined 
and smooth over Z^. 

In the sequel, any function f{xQ,Xc) will be understood as a function defined on Z (or on 
Zd by restriction) with f{s,9) = f{xo{s,9),Xc{s,9)) for {s,9) G Z. In particular, the function 
A'^ ()3.15p takes the expression 

1 3 

A-{s, 9) = e— ( ^ ^-PW\is,e) {r^-/y\ 0^ §). (3.17) 
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As stated in section [2l {A'^) is uniformly bounded in by 

\\{An-'\\oo < ^e^^ll'^ll--^ (rmn\^{e,e)\) ' < oo, (3.18) 

independently of oj. {A''^)~^ is also in L^^Zd), hence in every L^{Z£)) by interpolation. 
3.3.2 Radon Transform and Filtered-Backprojection operators 

In parallel geometry, for functions / and g respectively defined on M? and Z, we define the 
following Radon transform P and backprojection by: 

P[f]is,0):= [ f = [ fis9^ + t0)dt, is,e)GZ, (3.19) 
JL{e,s) JR 

pi[g]{x):=[ g{x-9^,e) de, X eM."^. (3.20) 
Note that P and P^ are adjoint to each other in the spaces L^(]R^) and L'^{Z) since 

g{s,9)P[f]{s,e) dsd9= [ P^[g]ix)fix) dx. 
In a 2D setting, inversion of the Radon transform can be done by means of the following formula 

This formula is however practically never used, because it requires differentiating the possibly 
noisy data as it tries to reconstruct / at arbitrarily small scales and this generates instabilities 
in reconstructions. Rather, we introduce a regularized version of (I3.2ip based on the following 
observation that 

PHwf, I P[f]] = W,lf = ^J^^ e-«<l /(O de, (3.22) 

where ★ and * respectively denote 2D and ID convolutions in the x and s variables and Wi, = 
is an approximation of identity {J- ^ denotes the inverse Fourier transform 
in 2D). Following [llj in two dimensions of space, we fix a function ^> first, which is a cutoff 
function in frequency supported inside [0, 1]. From we define the filter Wh by 
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where b plays the role of a bandwidth. Notice that the function Wb in (|3.23p naturally satisfies 
the identities 

Wb{u) = b^wiibu), and w'^\u) = b^+''w^^\bu), n > 0. (3.24) 
For g defined on Z, let us define the filtered-backprojection (FBP) operator by 

p-^^%] = P^[wbigl (3.25) 
and note the following straightforward estimate when g G L°^{Z), useful for section [37^ 

< 27Tb\\wi\\Li\\g\\oo. (3.26) 

Using convolution properties of the operators (j3.19p and (j3.20p . we can also express the FBP 
operator as follows [11] 

P-I'^ = — /-i'^ o ptt, where r^'^f := 

3.4 Proof of theorem EH 

We recall that 

R^''[k] = Tr~'^\Trn[k] + T^^k]] = V^P-''" [{An-\[TrR[k] + T^^^k]]] 
From theorem 13.21 we have 

WT^k] + T^+[k]\\L^i{dxr) < Cu-\ (3.28) 

where C depends on ||o"||(^2 and ||/^||c2- Reparameterizing + r2:!^[A;] into the cylinder 

variables {s,9) yields a same bound as (|3.28p . but in the space L°°{Z). Hence, using estimates 
(f3J8]) and I^Mh . we get 

< V^2nb\\wihi\\{An-^oo\\TrR[k]+T^4k]\\^ (3.29) 
Combining inequalities (j3.28p and ()3.29p concludes the proof. 

4 Proof of theorem 12.41 
4.1 Auxiliary lemma 

In the sequel, we use the following lemma twice, which is an exact first-order stationary phase 
expansion of an ID oscillatory integral with parameters. This particular case, where the second 



y ) \mf] 



(3.27) 
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derivative of the phase function is uniformly bounded away from zero along one direction, is 
useful for reducing the dimensionality of a given 01 while keeping under control the possible 
singularities (i.e. degenerate critical points) of the phase function. 
Let us consider the following integral for A G M", n > 0, 



/(A) 



J-'v^{^,^)f(^x,X) dx, 



(4.1) 



where / is assumed to be compactly supported in some set of the form [— A^;, A^;] xB, B compact 
in R", and ip is smooth on the support of /. Assume further the existence of < $m < < oo 
such that 



< A) I < $A/, (x, A) E supp /. 



(4.2) 



The lower bound in (|4.2p ensures that for every X £ B, there exists a unique X{X) such that 
dx(piX{X), X) = 0, since dx^{-,X) is strictly monotonous. By the implicit functions theorem, X 
is a smooth function of A. Condition (j4.2p allows us to make the phase globally quadratic in the 
X variable (which otherwise would require the Morse lemma in a neighborhood of the critical 
point, the size of which we could not control). We thus write 



99(2;, A) = S{X)+aK{X) 



'q{x,Xf 



where 



(4.3) 



S{X) := v{X{X),X), K{X) := |54^(X(A), A)|, and a = sign(5^,(^) 
is clearly constant on supp /. The lemma is stated as follows: 

Lemma 4.1. Let a one- dimensional parameter- dependent 01 be given by /(A) in (j4.ip . where 
if satisfies the condition (|4.2p . Then /(A) admits the following decomposition I = Iq + Ij., with 



/o(A) := 



V ^ ) K{X) 
3 



|/r(A)| < — mm 

OJ 



-y 



dx'- 



•,A) 



LH 



i=0 



57 
dx' 



(•,A) 



(4.4) 

XeB. (4.5) 



LH 



Proof. Starting from expression (|4.3p and using Taylor expansions, we have 

2 



rj = sign(3; — X{X)) 



(x-X(A)) 



K{X) 



{ip{x,X)-SiX)) 



K{X) Jo 



(l-t)|C.(/9|(te + (l-t)X(A)) dt 
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From this relation, we can show that the support of /, now a function of (ry,A), is included in 
[-Ho,Ho] X B, where 

Ho := ^y^. (4.6) 



The corresponding jacobian is given by 

dv _ 1 f^\d^,M{tx + {l-t)X{X)) dt 



dx a/2K(A) _ t)\d^^^^\(tx + (1 - t)X{\)) dt 

and is clearly well-defined everywhere by virtue of the bounds (|4.2p and smooth. Changing 
variable x — t- 77 in ()4.ip . we get 

/(A) = e*-^W / e*--^W^/(x(77,A),A)^(7?,A) dr?. (4.7) 

We now write a Taylor expansion of fdr^x with integral remainder in the 7? variable (note that 
dr^x = 1 at r/ = and that x(0, A) = ^(A)): 



/(x(7?,A),A)|^ = /(X(A),A) + 7?£|- 



/(x(.,A),A)|^(.,A) 



(tr/) dt. (4.8) 



In order to make the second term of the RHS of (|4.8p compactly supported, we write the 
following: let x is a smooth function of ry, supported in [—2Hq,2Hq\ and identically one on 
[—Hq^Hq], then xf = /• Multiplying ()4.8p by x(r/), we rewrite it as 

ox 

/(x(r7,A),A)— =/(X(A),A)+r?g(7?,A)+/i(r/,A), (r/,A)GMxi?, where (4.9) 



A) := xiv) 



d 



/(x(-,A),A)|^(-,A) 



(t?/) dt, (4.10) 



/o 5^ 

M^?,A):=/(X(A),A)(l-x(r/)). (4.11) 
Plugging ()4.9p into (|4.7p yields the decomposition / = /q + /r-i + /r25 where 

/o(A) := e-^W /(X(A), A) / e^'^^^'W^ ^77 = e^f ( —\ " e'-six) f{X{>)A) ^ ^4_^2) 

Jr V ^ / ^^(A) 



where we used the formula jj^e^'^'^'a dry = e^" i holds for cj = ±1 and /3 > 0. The 

remainder integrals are expressed as 



«2 g«I^S(A) /■ ^ „,,,„2 



I,i(A) :=e*'^^W / 7?e^'^"-^WV5(ry,A) dr7 = -a- — / e^'^'^W V (^, A) dr?, (4.13) 

/,2(A) :=/(X(A),A) / e-'^^W^(l-x(r?)) dr?, (4.14) 
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where we used the formula rie'^^'^^^^^~ = i^„\c(^x^ dr^^^^"^'"^^~ to integrate by parts in /ri(A). 
We integrate by parts once again in (j4.13p and obtain this time 

IgiX) = a- — ^ / F^(ryv4Z^)-4(r?, A) dry, (4.15) 

■ (2 

where the function F„{u) := dt is related to the Fresnel integral and is known to be 

uniformly bounded in M. From the two expressions of Iri (A) in (j4.15p and in the last right-hand 

X 

side of (|4.13p . the expression of g in (j4.10p and the fact that — — is bounded for n < 4, we 

orf^ 

obtain that Iri satisfies a bound of the form (j4.5p . Moreover, the remaining integral Ir2 (|4.14p is 
rapidly decaying because 1 — x(^) is supported away from = 0, so it can be made comparable 
to (if not smaller than) I^i- Indeed, we have for every p > 



p 





' d 1" 




drj rj 



p 



[l-xMdii, (4.16) 



where we have used p integrations by parts. From equation (|4.16p with p = 1, we obtain a 
bound of the form 

|/r2(A)| < C|/,i(A)|, AGS, 

where the constant depends only on the features of the function x- As a conclusion, setting 
Ir := Iri + Ir2, it is clear that Ir satisfies bound (|4.5p . The proof is complete. □ 

4.2 Main outline and proof 

In order to prove theorem 12.41 we must bound the operators R!^' in ([211]) and R^' in (f2l2]) 
and make them contraction operators in some spaces to be defined. This is achieved by doing 
stationary-phase analyses of the compositions of operators instead of bounding each operator 
separately as was done when proving theorem 12.21 

The proposition below bounds the error operator i?"^'^ in (|2.1ip and shows how close the 

operator Tf ' comes to providing [pk]^. 

Proposition 4.2. The operator R'^"^[k] is bounded in C{L?{Br-2D)) and C{L°°{Br-2D)) and 
satisfies the following bounds for a; > 1 

\\Rf\\ciLHB^^,,))<^{b + b') (4.17) 
\\RT''\\ciL^iB^.,o))<-ib' + b'). (4.18) 
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Let us now introduce the operator L (first introduced in [T]) defined as: 

Lu{x) = [ }ME^yl dy. (4.19) 
Jx \x-y\ 

It is proved in [7] for instance that L is weh-defined in C(U'), p G [l,oo]. Its operator norm in 
each of these spaces, denoted by ||-/^||p, satisfies the estimate 

E{x,y) 



imip ^ 11-^11 '■= sup / dy < 27rA. 

x&xJx \x-y\ 

The next proposition controls the second error operator i?2'^ (j2.12p . provided that the scat- 
tering coefficient is uniformly bounded by (||(/'||oo||-^||oo)~^- 

Proposition 4.3. Let < Kq < (||0||oo ||^||oo)~^ and k,k e i3/<o(L°° (5^-21))) ■ Then, there 
exists a constant C2 that depends on \\a\\c2, Kq and A such that, we have 

Rf[k] - Rf[k] < C2-^ log (^) \\k - fclloo. (4.20) 



00 



UJ2 



Remark 4.4. The boundedness conditions for k, k in proposition \4-3\ are used to ensure that 
the scattering series converges. They are sufficient but not necessary. A less stringent condition 
would be to impose that the operators L^k and L^j^ have operator norms in C{U\X x S"*^)) 
bounded by 1, where we have defined 

Ta,kf{x,v)= [ k{x')j^ — ^—^(j){x — x',v)f{x,x — x')dx, (x,u)sXxS^. 



X 

For the sake of clarity, however, we do not consider such a proof here. 

The proofs of propositions 14.21 and 14.31 above are given consecutively in the next two subsec- 
tions. We now prove theorem 12.41 

Proof of theorem \2.4\ Let us fix a reconstruction bandwidth b and < Kq < (||</'||oo||-^^||oo)~^- 
Applying the inversion to the data yields the equality 



Tf ^'"D" = [kp]i, + R^'^ikp] + R^'^ikp] := [kp]b + R''''[kp]. (4.21) 

The introduction of p into the argument of the remainder operators may change the constants 
Ci, C2 of ()4.18p and ()4.20p but causes no trouble since p is bounded from above and away from 
zero on Bj—2D- The bounds (|4.18p and (|4.20p decrease to zero as c^j — )• 00 so that there exists 
cjo such that for uj > ujq, we have 



^(62 + 6^) + C72-^logm=ci<l. 
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Thus, when uj > ujq, the operator R^'^ = R'^' + R2' is a contraction on Bko{L°° {Br~2D))- We 

must now show that the mapping k 1— )■ T'f ' D'^ — R^'^(k) is a contraction in Bki{L°° {Br-2D)) 
for some < Ki < Kq. In order to do so, we write, using (j4.2ip . 



^'''ti^ - R'^'''[k] = [kp]b + R'^^^kp] - R'^^^k]. 



For a given x G Br^2D, we have 



[kp]bix)\ 



Wbiy - x)[kp]{y) dy 



< llWbllill^plloo, where 



Wbix) = j^^ e*"-«4> (^^^ di = thus WW^Wl^ = \\W^\\l^, 

and so we obtain the bound 

||Tf ~''V - R'^\k)\U < {\\Wi\\li + ci)||A:p|U + cillfelloo. 

Therefore, setting Ki := assuming that ||A;/9||oo < Ki, this ensures that ||fe||oo < 

Kq imphes \\Tf ' 2)'^ — R'^'^(k)\\oo < -K"o- Hence, whenever ||A:/9||oo < -f^ii the iterative scheme 
(j2.13p converges to a hmit point k* G BKQ{L°°{Br-2D)) tbat satisfies 

Tp^^'^D"^ = k* + (4.22) 

Subtracting ([02]) from (fOT]) yields 

r - [Mb = ^'^'iM - ^'^'1^1- 

Thus, we have 

W - [Mblloo < - i?"''[r]||oo < CiWkp - r lloo 

< ci||/cp - [/c/o]b||oo + ci\\k* - [kp\b\\oo, hence 

- [/C/9]b||oo < — \\kp - [kp]h\\oo- 
i — Cl 

Theorem 12.41 is proved. □ 
4.3 Proof of proposition 14.21 

The operator ii"^'^ is hnear in /c, and we now exhibit its Schwartz kernel. We recall that 
T^[k]{s,9)= I e-^"(l"-"°l+l"-"-l V(^^^o, ^7^x)E{xo, x, x,) '"^ ~ " ' ~ f ' [pk] (x) dx. 

Jm? \x — Xq\\x — Xc\p[X) 
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We now apply to it the inverse operator ()2.9p . using the expression p.27p of the FBP operator. 
We first have 

{A'^)-^xTi[k]{s,e) = J— [ e*"'^^(''^'^')c(s,6',x)[pA;](x) dx, where 
V 27r J^2 

(/7i(s, 9, x) := \xc — xo\ — \x — xo\ — \x — Xc\, (4.23) 



c{s,9,x) := xis)xi\x\) 



(^(ar^o,xr--^) E{xo,x,Xc) laT^o • Uxo\\x^-~Xc ■ u^J \/2r^ 

3 
)4 



ip{0,0) E{xo,Xc) \x - Xo\\x - Xc\p{x) (^2_g2): 

(4.24) 

Applying now a backprojection and exchanging the integral signs, we write 
V^P^iiA^r^T^my) = [aiy,C,u)T[pk]] (y), where 

aiy, uj):=J—[ [ e^(^-^>«e*"^i(j'-^^'^'^)c(y • 6^, 6, x) dO dx. (4.25) 

V 27r J^2 Jgi 

Finally applying the operator I^^'^ in (|3.27p . the operator T'^ ' can be expressed as 

TT"'" V[fe](2/) = -^if^/''^^ (y) \iHy.^.u^mpm) d^, (4-26) 

with a defined in (I4.25p . 

Proof of proposition \4-S\ We now analyze the symbol a{y^ ^, u) in (I4.25P by studying the behav- 
ior of the phase function ipi ()4.23p in terms of x. For fixed y € Br~2D and 9 £ S^, we have 
VxVi = whenever 

X — Xq + X — Xc = 0, 

which happens exactly when x € [xo,Xc]. It is thus convenient to use the change of variable 
X = y + u9 + v9^ for (w, v) G M^, and the symbol a becomes 

a{y,tc^) = \h^ [ [ e'<-'l{y,^,u,9)dud9, where (4.27) 

V 27r Jgi 

I{y, e, u,9)= [ e''-MyO\o,u,v)f^^^^ ^428) 

fx{y,u,y,9) := e'<-'\{y ■ 9^ ,9,u,v). (4.29) 
We are now in the setting of section 14.11 where the amplitude function fi is supported in 
supp /i C ^{v,u,y,9) G X Br^n x s.t. x{u,v,y,9) G Br-D and \y ■ 9-^\ < r - 

(4.30) 

C [-A, A] X [-A, A] X Br^D X (4.31) 
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On the other hand, the phase function, expressed as 



\/v'^ + \Px — xop — \/v^+\Px 
is proved in appendix El to satisfy 

-OO < -^>l,A/ < < -^l,m < 0, with 



D 



{Px ■.= y + u9) 



4^2 



(4.32) 



uniformly on supp /i. Moreover, since (/?! is clearly even in the unique v that satisfies d^^i = 
is = regardless of all other variables. Hence, section HTTl gives us the expression of a rigorous 
decomposition for the integral / in (j4.28p . Note that in this case, we have 

Si{u,y,6) := ipi{y • 6^,6,u,0) = 0, Ki{u,y,6) := d'^^ipi\y=o = -dop{y + uOf . 

fi{0,u,y,9) = y^p{y + u§). 

Thus the integral / in (j4.28p decomposes into Iq + Ir, where, following formula (j4.12p . 



h{u,y,0) = e 



<jj J " Ki{u,y 
The remainder term Ir satisfies the following bound, following (j4.5p : 



\I,.(u,y,9)\<^Yl 



i02 



1=0 



97: 
dv 



U-,u,v,9) 



(4.33) 



From expression (I4.29p . it is clear that for < i < 3, 



d'fi 



which in turn implies the bound 



j=0 



C 



\Ir{u. 9,y)\<—{l + leh, {u, 9) G supp c. 
a;2 



(4.34) 



(4.35) 



Plugging the decomposition I = Iq + Ij. into expression ()4.27p yields the following decomposition 
a = oo + where 



ao(y, ^)-=\l^ I / e™«'^^o(^x, y, 9) du d9 = j / e*"«-^ du d9 = ' 2 



27r 



§1 JR 



§1 



47r 



R2 \A \i\ 



ar{y,(.,u}) := \j — 



z'<-hr{u,y,9) du d9. 
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By virtue of bound (|4.35p , we get the following bound on : 

\ar{y,tio)\<-{l + \^f). (4.36) 

Plugging a = ao + Qr into expression (j4.26p shows that Ti[k]{y) is the sum of two terms, 
the first one being 



L (t) I*o(^.«.-)^M(«) « = X, -"'i if) 

and the second one being the error operator 

Taking norms, we obtain the first bound 

/ \Rf[k]{y)\Uy = Ci [ ^(^f)'\^\'\ar{y,^,u)\'\T[pkm\' 

-^'^L^ (y) (1^1 + i^i'fmpkmf d^ 



npm) dc = [pkUy), 



<c,^{b+b^f [ \j^[pkm\'d^, 



where we used that < <l> < 1 and that on the support of <l> (4^), we have |^| < b. Using 



^ b 

2 



Parseval on the latter right-hand side, this yields the C{L ) continuity inequality 

\\R^Amv^<Vc'i-{h+h')\\pkU2, 



and hence the bound in (|4.17p since p is uniformly bounded on Br^D- 
In order to get an C{L°°) bound, we have the following estimate 

<T^(''= + i>')MlL»ll*lln, 

where we used the following estimate based on the Cauchy-Schwarz inequality: 

l^rd \Fm)\ di < \\gh.\mL.b'>+\ q > 0, 

and hence the bound (|4.18p . This completes the proof. □ 
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4.4 Proof of proposition 14.31 

For a, /3 formal subscripts, denote Va^p := xp — Xa- For m > 2, denoting = {xi, . . . , Xm) and 
= dxi . . . dxm., the multiple scattering of order m is given by 

J '^J — 1| \Xm Xc\ 

(4.37) 

We now multiply pointwise by {A^)^^{s, 0)x{s) and apply the FBP operator using expression 
(|3.25|) . The operator ' can be written in the compact way 

f^~^'''Ti^[k]{y) = V^ [ /;;;[fc](x„)/3"(y,a;i,a;2,x^-i,x^) d^m, (4.38) 

where we have defined 

fuJHTf ^ r f \uf .E'^{Xm-l,Xm)"fr^ Hxj)E'^iXj-l,Xj)<l){Vj_ij,Vjj+i) 
/mW(Xm) = k{xi)k{Xm.)-. p J_J_ ^ '■'^-^ , 

(with the last product replaced by 1 when m = 2), as well as 

fS'^iy, xi,X2,Xm-i,x,n) = [ f e*'^'^2(^'^'"i'"'")u;fe(2/ • 6^ - s)a"(s, 6, xi,X2,x,n-i,Xm)xis) ds de, 

(4.39) 

a'^{s,9,Xi,X2,Xm-l,Xm) ■= {s,i 



X Wxo -VoMWxc ■ Vm,c\(l){vO,l,Vl,2)4>i'"m~l,m,Vm,c), (4.40) 
ip2{s,0,Xi,Xm) := 1X0(5,6*) - Xc{s,6)\ - 1X0(5,6*) - Xi| - |Xc(s,6') - Xm\- (4.41) 

Since we do not assume any regularity on k, we cannot use stationary phase with respect to 
a variable on which k depends. This is why stationary phase arguments will be used to bound 

(3'^ but not the /^[/cj's. This requires estimates on /3'^ in order to bound the operator ' 
using Holder inequalities. We will use the following: 

Proposition 4.5. The function admits the decomposition (3^ = /3^ + /32 , where the functions 
I3f and satisfy the following bounds for every {y, xi, X2, Xm~i,Xm) G {Br-o)^ such that 

Xl 7^ X-m- 

— 1 

|/3^(y,xi,X2,Xm-i,Xm)| < I3f{xi,xm) ■■= Cp,i . , (4.42) 

W yj\xi - Xm\ 

— fet 1 

\l32{y,Xi,X2.,Xm-l,Xm)\ < /3^(xi,Xm) := C«.2 -. r. (4.43) 

iO Xl - Xm\ 
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The proof is rather technical and is postponed to section 14.51 We are now in a position to 



bound the operator R'^''^ (j2.12p . 

Proof of proposition We start from equahty (j4.38p and obtain estimates for Let k,k £ 
L°^{Br~2D)- Then using the fohowing property that holds for real numbers ai, . . . , Um, bi, . . . ,bm 



ai - ■ ■ am - bi - ■ - bm = '^bi - ■ ■ bj^iiaj - bj)aj+i 



(4.44) 



with aj = k{xj), bj = k{xj), j = 1 . . . m, and the fact that max(||fe||oo, ||fc||oo) < Kq, we have the 
following estimate 



|/-[A:](x)-/-[A;](x)| <mi^o" 



m~l II I ||m— 2 



n 

i=2 



E{xj_i,Xj) 



\k-k\ 



(4.45) 



Using estimates (|4.42p and (|4.43p together with (|4.45p . we have that for k,k e BKo{L°°{Br^2D)) 
and m > 2, 



m-2 I 



k 



)[Crrt,l + Cm,2], 



(4.46) 



where Cm,a ■= /3'^{xi,Xm)Pm{xi,Xm) dxi dxm, a = 1,2, (4.47) 

m 



Pm \ Xl , Xui j '.- 



Y^ E{XJ.^,X^, dx 



^"^^ j'=2 ~ 



m — 1 



L 



m-2 



E{xi,-) 



\Xi 



In order to bound the Cm,aS, we need estimates on pm- We have the following estimates 
(obtained by bounding the attenuation terms by 1 and performing polar changes of variable 
inside the integrals): 

P2{xi,X2) < \x2 - xi|~\a.e. (xi,X2), 
P3{xi,X3) < Csl log(|x3 - 3;i|)|,a.e. (xi,X3), 
Pm{xi,Xm) < C4\\L\\^~^,a.e. {xi,Xm), m > 4. 

Hence by virtue of the bounds ()4.42p and (|4.43p . the only term that requires attention is C2,2- 
The integrals defining all the other Cm,a are well-defined. Addressing 6*2,2 1 we change variables 
Xm = xi + pa and split the integral in p into two intervals [0, e] and [e. A] and use different 
bounds for /Jg as follows 



C2,2 < Cp,2 



< Cp,22tt\X\ 



Ix J\ 


r r 
51 JO 








e + 







dp da dxi + 



1 

6\ 2 



dp 
P 



da dx 



6\ 2 



UJ 



(log A - loge) 



<Ci,^logf, 

' oj b 
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if we choose e = (^) ^ • For the other terms Cm, a, we obtain bounds of the form 

C2,i < C2 1 — , C'3 q,<C3 — , a = 1,2 

' 00 ' 00 



n ^ n " \\T ||m-4 
'^m.n ^ '-'4 ll-'^lloo ' 
00 



m,a ^ ^4 — 11-1^11^ , a = l,2, m>4. 



Hence for m = 2 we have 



while for m = 3 we have 



\k-k\ 



<2,KlU\\[C',^,+C', 



3,2J 



(4.48) 



(4.49) 



Summing inequahties in (j4.46p for m > 4, we get 

00 

Y,\ff^T^[k]-Tf^T;:i^> 



m=4 



< 



,3 00 

2C'4^^^II^ - ^lloo V m(i^o||</'|| Halloo) 



m-2 



< 2Ca 



m=4 
1 63 



|L||2,(l-Ko||'^||||i|| 
where we used the fohowing inequahty for |x| < 1: 



\k-k\ 



(4.50) 



E 



1 



(1-X)2' 

Finally using the fact that 

00 

Rf[k] - Rf[k] < ll^r"'' - T^'^'^TZik] 

m=2 

we obtain the desired result using bounds (j4.48p . (j4.49p and (|4.50p . □ 
4.5 Proof of proposition 14.51 

The function (3^ is an 01 in {s,9) parameterized by the points {y,xi,X2,Xm~i,Xm,) G (-Br-2Z))^5 
which we now analyze using stationary phase techniques. The only parameters that matter for 
the analysis are those the phase function ip2 in ()4.4ip depends upon, that is xi and Xm, and all 
estimates are uniform in (y, X2,Xm-i)- Critical points of 1^2 become degenerate as l^i — Xm\ — ^ 
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and so a careful analysis is necessary to control this degenerate behavior and capture the leading 
term as — t- oo. 

Let us analyze the phase function ip2 in (j4.4ip first and in particular which couples (s, 0) 
satisfy {ds'^2-,deip2) = (0,0). Using identities from section [Bl the gradient of ip2 in the {s,9) 
variables can be written under the form: 

r 

dsV2 = —F==={k- fc), and deip2 = r{fo + fc), where 
V T — 

/o := {xi^o - 0) ■ dsXQ, fc := (x^^m - 0) ■ dgXc- (4.51) 

Now, since \J r"^ — is never zero on 2/5, the system of equations {ds^2-,de^2) = (0,0) can be 
recast as 



- s'^ds'^2 ± d0ip2 = 0, 

which is equivalent to fo = fc = 0. This in turn implies xi — xq = Xc — Xm = G- Hence the 
critical points of ip2 are exactly those for which the points xo,xi,Xm,Xc are aligned, that is: 

if Xl / Xm, two distinct points (s, 6) G {{xi ■ 6im , Oim), (xi • 0im + , dim + tt)}, where we 
defined 

:= arg(xm - Xl). (4.52) 
At these points, we have (/?2 = \xi — Xm\ and —\xi — Xm\, respectively, 
if Xl = Xm, a curve in Zd of equation {s = xi ■ 9^}, along which (p2 = 0. 

The second case describes the caustic set in the space of parameters, i.e. the parameter values 
for which the stationary points are degenerate. Since (p2 is constant along this curve, it is clear 
that these stationary points are degenerate, see [121 [13]. This degenerescence can be observed 
by noticing that the determinant of the Hessian at the two critical points in the first case is 
proportional to \xi —Xm\ and hence goes to zero as \xi —xi \ — )• 0. The degenerate critical points 
are not of type Ak (as defined in [T]) for any k > 2, since the phase is locally constant along 
one direction as one passes through a degenerate critical point. However, because at any of 
these critical points, rank H^^ > 1, the approach of |13[ Section 3] which we are following here 
remains efficient in order to control the degenerescence near the caustic set. 

Following this approach, we first need to find a direction along which the second derivative 
of is uniformly bounded away from zero. In our case, it is proved in appendix IB. 21 that for 
any {s,9) £ Zd, 

-oo < -<I>2,M <^^is,9) < -<^2,m <0, with $2,m = ^ and $2,M = -^. (4.53) 
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As a result, for every 9 € S^, the function dsip2{-,9) is a strictly decreasing function. More- 
over, it is easy to check that dsip2{xi • 6^,9) and ds^2{xm • 0^,0) have opposite signs. This 
ensures that for every 9 G the equation ds^2{s,9) = has a unique solution called cr{9), 
which satisfies the following bounds 

min(xi • 9^,xr,^ ■ 9^) < a{9) < max(xi • 9^,x„i ■ 9^). (4.54) 

Except for the simple case where xi ■ 9^ = ■ 9^ = cr{9) (i.e. xo,xi, Xm-,Xc are aligned), it 
appears that the computation of (7(9), which requires solving the equation /o = fc (with /o, fc 
defined in (|4.5ip ) leads to quite complicated expressions. It is however important to note that, 
because cr{9) satisfies the following ODE (derived from the implicit functions theorem) 

da ^ d%^{a{9),9) 

d9 dl^{a{9),9y ^ ■ > 

and because ^2 is smooth with respect to all of its arguments, then u{9) is a function of 9. 
For the sequel, we now define the functions 

S2{e) := M^{O),0), K2{9) := \dl^2{cT{9),9)\ = -dlM<0).0)- (4-56) 

Let us now assume x\ 7^ x^. First notice that, using the chain rule as well as relation (j4.55p . 
the first two derivatives of ^2 (|4.56p are given by 



S'^{9) = dem{<0).O), and S'i{9) 



detH^,{a{9),9) 



K2{9) 

The critical points of 52 (i.e. the zeros of ^2) are 9\m and + vr, at which we have 

'52(6'im) = -| 1 — I r, and 5 ((9im + vr) = — — -. (4.57) 



Proof of proposition Let us denote by 

(3'^{y,xi,X2,Xn,.i,x^)= [ I2{9)d9, where l2{9) = [ e'^'^^ f{s,9)x{s) ds, (4.58) 
f{s, 9) := a'^{s, 9, xi,X2, Xm-i,Xm)wb{y ■ 9^ - s). (4.59) 

The dependency on parameters is made implicit to improve readability. Using the identities 
(|3.24p satisfied by Wh and the fact that a"^ is smooth, we have that for every < i < 2, there 
exists a constant Cf^i such that 



<Cf^ib{l + b'), 9€§\ (4.60) 
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Let us focus on the 01 /2(^) in (|4.58p first. Again, since we have the property (j4.53p . we can 
use the results of section ITTl We can decompose I2 = hfl + l2,r, where foUowing ()4.12p . 



(27T_Y fia{e),e) 



K2{e) 



c 



j=0 



57. 



< ^6(1 + 62) 



LH 



UJ 



(4.61) 
(4.62) 



Integrating now hfii^) with respect to 6, we get an oscillatory integral J2,o := /gi hfii^) dO 
with phase function 82(0). We recall that the critical set of 5*2 is ©c = {Oim, &im + tt}. We now 
claim the following 

Lemma 4.6. There exist constants | > > 0, Ci > and C2 > independent of {xi,Xm) 
such that the following inequalities hold: 



\S'2{6)\ > Ci\xi - Xm.\ for min 16* - 8c\ > ^, 
|5'9(6')| > C2IX1 — Xm\ for min \6 — OA < Sq. 



(4.63) 
(4.64) 



Estimate (j4.63p will help us bound the integral J2,o away from the critical points, while 
estimate (I4.63P will help us bound the contribution to J2,o in a- neighborhood of the critical 
points. Lemma 14.61 is relegated to appendix IB. 4[ Let 60, Ci and C2 be given by lemma HTHl Let 
^1(6*) be a smooth function such that pi{9) = if |0| > 60 and pi{0) = 1 if \0\ < ^. We now 
write the following partition of unity 

l=p,{6) + p{e), where p,(0) := pi(e - 0,), 

See 0c 

thus p{6) = if miug^g©^ \d — Oc\ ^ Using this partition of unity, the integral J2,o can be 
written as 



-^2,0 = 6 *4 ( — 1 [J2,o,i + -^2,0,2] 



where 



J- 



2,0,1 



i^2(^)2 ' ' V§1 ^2(^)2 



The integrand in ^2,0,1 vanishes for min^^ge^ \9 — 6c\ < 5o and so we can use estimate ()4.64p 
together with [H Corollary 2] to bound ^2,0,1^ 



\J'. 



2,0,1 



< 



12 



y/ujC2\xi - Xr, 



d_ 

de 



—Pc{0) 

K2i9)-2 
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— ^201 

ii(si) ^/u;\xi-xm\ 



(4.65) 
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Treating J2,o,2 now, a first bound can be obtained easily 



\J: 



2,0,2 



< 



K2{ey- 



Li(§i) 



< ^202^ . (4.66) 

Then, using the fact that the integrand in J2,o,2 vanishes for 9 such that ming^gQ^ \6 — Oc\ > ^, 
we can use the bound (|4.63p to integrate by parts: 



•^2,0,2 

whence the bound 



1 


1 - 











f{<o))~p{e) 



2,0,2 



< 



UjCl\xi - Xmp 



K2{e)-2S'^[e) 



de = — [ e-^^w^ 

./si dO 



S'2{0) 



.d_ 

de 



[K2{9)-2S',{e)\ 



< Co 



de. 



Taking the geometric average of bounds (j4.66p and (14.67P , we get 

bl 



\J2,0,2\ < (C'202C'202 



(4.67) 



(4.68) 



UJ2\Xi — Xr, 



To conclude, if we define /5f := e''^ (^) ^ J2_o,i + /gi hA^) ^9 and := e"*? (^) ^ J2,o,2, we 
obtain the desired decomposition. Estimate ()4.42p is obtained by using inequalities (j4.65p and 
setting := \/27rC2oi + 27rCj,, while estimate ()4.43p is obtained by using inequalities (14.660 
and ()4.68p and setting C/3,2 := \/27r max(C202, (C'202C202)^)- 



□ 



A Formulas for proposition 14.21 

From equality 

ipi = \Xc — Xq\ — y/v"^ + \Px — XoP — + \Px — Xcp, 

we get that 



\Px — Xq\^ ^ \Px — Xc 
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\x — Xo\ 



\X — Xc\ 



Thus we have 

\Px — XqP \Px — 



\x — Xo\ 



+ 



1 



\X — Xr 



3 >^{\Px-Xo\^ + \Px-Xc\'^)> 

\xr) — xA'^ r"^ — D 
> L-^ _£!_ = _ > 



{\Px — xq\ + \Px — Xc\y 
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and hence the lower bound on |(9^„(/?i|. The upper bound is given by 

IPX — XnP \Px — xJ'^ 1 12 



\x — xo\^ |x — Xcp |x — xqI Ix — D 

B Formulas for proposition 14.51 

B.l Identities involving the boundary points xq, Xc 

Concerning the boundary points in the Radon geometry, we have 



OS _ §2 \ V r j as yjr^ _ ^2 \ v r 

The second partial derivatives satisfy the relations: 

9^X0 _ d'^XQ _ 1"'^ a d'^XQ _ 1 

d'^Xc d'^Xc —r"^ X d'^Xc —1 

B.2 Properties of 

Derivatives of the phase function: we have 

f2 = \xo - Xc\ - \xo - Xi\ - \Xm - Xc\, 
ds(p2 = , ^ ^ „ (2- +XO-XI- dsXo + Xc-Xm- OgX, 

r 



^ {(xx - Xq-O) ■ dsXo - (Xc - Xm-0) ■ OgXc^ 



'r^ — s 

= r ( xV^i • deXQ + xT^m ■ dgx. 



r (^{xi - xq- 0) ■ dgXo + (xc -Xm-0) ■ dgXc^ 



In order to compute the second derivatives, we use the following identity, where x depends on 
(o, b) and y is a constant vector: 

, , — d^x 1 / dx dx ( dx\ f - — - dx 



dadb^^ y ^ y Qadb^ \x — y\ \da db \ ^ da) \ ^ db 
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We obtain 



dsde - s2 

5V2 



(r^ — s^^ 2 '- 
1 



a^cT-^i • a;o + xT^^^ • Xc + r2(/o - fc) 



d9^ 

where we have defined 



Xq - Xi ■ Xq + Xc - Xm ■ Xc - r (/q + fc), 



1 - (x^r^i • Ss^o)^ , „ 1 - {x'T^rn ■ dgXcf 

Jo ■■= , , and fc := . . . 

\Xi Xq\ \Xc Xrn\ 

Since xq (resp. Xc) is orthogonal to dgX^ (resp. dgXc), we can rewrite /o and fc as 

1 (X(P^1 • xo)^ , „ 1 (x^^m • Xc)^ 
Jo := — \ \ and fc := — | ■. . 

|xi — xqI \Xc — Xm\ 

In order to show estimate (I4.53p . the fact that xi and Xm belong to Br^2D yields the bounds 



, -—^ ^ ^ j2D(2r - 2D) 2D 

min(xo - xi • Xo, Xc - Xm ■ Xc) > > \ . 

r V r 

Together with the obvious fact that 2 + (xq — xi — Xc — x^) • ^ > 0, we obtain that 



9^2 

Qg2 — j^2 _ g2 y _ 2;q| \Xc — Xm\ J T f 



< 



1 1 

7 + 



2D 2D 
— <- — 



and hence the upper bound in estimate (|4.53p holds. The lower bound is straightforward if we 
use the following inequalities — > rD, |xi — xo| > -D and \xm — xd > D. 

B.3 Determinant of if^^ at a point {a{6),6) 

This appendix is useful in order to compute 5*2 at the critical points. The condition ds(p2{s, 0) = 
that defines cr{0) can be rewritten as 



s(2 + (x'JP^i - xT'^rn) ' &) + \/ r"^ - s'^{x^^^i + X^^m) • (9"^ = 0, 

which in turn allows us to obtain the relation 



2 + (xo - Xl - Xc - Xji 



2/^2 _ ,2 



- - Xo - Xl • Xo - Xc - Xm • Xc 



(B.l) 
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Defining 



go := xq-xi-xq- r^/o = xq- xi- xq\1 



dc ■ — ' ^ fc — ' I 1 



FO ~ Xi\ I pi — a^ol 

ry -y» , \ I rf ry \ I ry ry , ry \ 

I rf rf I I I /v" fji I 



and using relation (jB.ip . we can write the second partials of (p in the following compact way: 
-^ = ^T^ 5o + 5c-2VH-s' , ^^=50+ 5c, = ^==(5c - 50), 



thus the expression of the determinant of the hessian at a point {(j{9), 9) takes the form: 

= (2»* - v^to + *)) = - i) - i) - 1. 

This expression is zero whenever x\ — x^. 
B.4 Proof of lemma [4.61 

We first claim that there is a constant such C3 such that, 

\S'^m<C2,\xi-Xm\, 9e8\ (B.2) 

Indeed, the function {9, xi, p, a) ^ S'^'{9, xi, xi+pa) is smooth with respect to all its arguments, 
and such that limp_^o 'S'2" ~ 0, so we have that 

S2(9,xi,xi + pa) = p —-^{9,xi,xi + upa) du, 
Jo dp 

dS'" 

and since — - — is continuous on a compact set, it is uniformly bounded by some constant C3, 
op 

hence the estimate ()B.2p . Using a Taylor expansion of 5*2 (0) about 9(. E 0^, we thus have: 

\S2iO) - SUOc)\ < \i0-9c)\ C\S'^{u9 + {l-u)9c)\ du<\9-9c\C:i\xi-Xm\. 

Jo 

Hence for \9 — 9c\ < I'^z^^'^^l -^^e have that 

\S'm\ > hsm\ = ^^^^ rki - > - x^\. 

2 2\xi - xo\ + \xc - Xm\ or 
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Hence, if we define 60 
Co := 1 ^ 



1 V2rD 
2C3 4r 



< 2C^min0^6ec 



\suec)\ 

\xi-Xm.\ 



, condition ()4.64p is fulfilled with 



We now focus on proving inequality (|4.63p . Let us define the function 

G{s, 9) = (r2 - s^)id,ip2f + {deip2f = + fl 



and notice that |5'2(0)| = G{a{9),6)'2 . Let us fix 9 such that mingi^ge^ \9 — 9c\ > By virtue 
of bounds (j4.54p . there exists A € [0, 1] such that cr{9) = (Axi + (1 — X)xm) ■ so that 



Xl-XQ= [(1 - \){xi - Xm) 



+ [\/r2 - 0-2 + xi • 9]9, 



Xr Xr 



[A(xi - Xm) ■ 9^]9^ + [\/r2 - a2 - x„, • 9]9. 



Treating /o first, we split cases: if (xi — xq) ■ 9 < 0, then necessarily |cr(^)| > D, and cr{9) and 
(xi — Xq) ■ 9^ have opposite signs. Hence, 



l/ol = (1 - xi - xo • 9) h \xi - xo 

r 



a' D 



Now if (xi — Xq) ■ 9 > 0, we have 
1 



l/ol 



where 



\Xi - Xo\ 



{Xi - Xq - \xi - Xo\9) ■ dsXo 



{I - X)\{x^ - Xm) ■ 9^\ 



\Xl - Xo\ 



a 

1- — -7- 



7 := 



Hence, if cr < 0, then 



and if o" > 0, 



- Xq\ - (Xi - Xq) 
(Xi - Xo) • 9-^ 



(Xi - Xo) 



|X1 - Xol + (Xi - Xq) 



< 1. 



> 



l-Zo -7- 



+ 



Finally, we have that 

|(xi - Xm) 



xi - Xmll sm^yim 



, a 1 
+ (l-7)->2- 



> |xi - Xm| sin—. 
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In summary, denoting C := min sin we obtain the bound 



l/ol > min ( ^, (1 - A)|xi - Xm\C 



Likewise, we obtain the bound 



which imphes that 



|/c| > min ( X\xi - Xm\C 



= V/o+/c > min (7, - ^-1) 



where we used the inequahty + (1 — A)^ > ^ for every A G [0, 1]. Lemma HT6l is proved. 

Remark B.l. Although we do not know the function S2 explicitly, we know that we can not 
expect a uniform bound of the form. 

\S'^{0)\>C\xi-Xm\, 9eS\ 

(which in turn would give us a less singular bound of the form |/3'^| < C\xi — Xm\^^ by virtue 
of m Corollary 2]). Indeed, since S2 has two zeros and is periodic, its derivative has at least 
two zeros. In addition, away from the critical points, as \xi — x^l — 0, ^2 cannot decay to zero 
slower than \xi — Xm\- This is because we can obtain an upper bound for jSg] similar to ()B.2p 
for 1 5*2" I . This in turn guarantees that \S2\ — )• at least as fast as \xi — Xm\ does. In that sense, 
lemma \4-6\ cannot be improved. 
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